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Abstract 

The temperature distribution in the nuclear fuel rods for high burn- 
up is studied. We use the numerical and analytical approaches. It is 
shown that the time taken to have the stationary thermal regime of 
nuclear fuel rod is less than one minute. We can make the inference 
that the behavior of the nuclear fuel rod can be considered as a sta- 
tionary task. Exact solutions of the temperature distribution in the 
fuel rods in the stationary case are found. Thermal regimes of high 
burn- up the nuclear fuel rods are analyzed. 

1 Introduction 

Temperature distribution in the fuel rod of nuclear reactor is one of the 
most important factors that controls the behavior of fission products in the 
pellets, the diffusion and vaporization properties and so on pQ. This depen- 
dence has been intensively studied for prolonged lifetime of existing reactors. 
There are many papers in which the authors studied the different aspects 
of this problem [2H7]. Yapici et al. [2] investigated the maximum tempera- 
tures in centerline of the fuel rod for different clad outer surface tempera- 
tures, melting points of the fuel materials, temporal heat generation, tem- 
perature distribution in the nuclear fuel rod and temporal variation of the 
neutronic data during rejuvenation periods. In [3J Kim et al. developed 
the one-dimensional heat conduction model to determine the temperatures 
distribution from the fuel center to cladding surface in the radial direction. 
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Pontedeiro et al. in [I] presented an improved lumped — differential formu- 
lation for one — dimensional transient heat conduction in a heat generating 
cylinder with temperature — dependent thermo — physical properties typical 
of high burn-up nuclear fuel rods. Analytical model for the determination of 
the temperature distribution in cylindric heater components with character- 
istics of nuclear fuel rods is given by Fortini et al. in [5]. Espinosa — Paredas 
et al. in paper [7] explored the applicability of a fuel rod mathematical model 
based on Non — Fourier transient heat conduction as constitutive law for the 
Light Water Reactors transient analysis. 

The fuel behavior is affected by the temperature distribution in the fuel 
that is related to change in the fuel microstructure with irradiation [18]. One 
significant change in the fuel microstructure is the formation of a porous 
rim in the periphery of the high burn- up fuel [8"HT7j. High burn-up nuclear 
fuel rods have been intensively studied last years [I]. Many papers were 
published studying of the rim formation mechanism [T8T[20] . It was shown 
that the rim structure is formed through recrystallization and coarsened pore 
formation. The subdivided grains with high angle grain boundaries are the 
nuclei of recrystallization, and then the coarsened pores are formed by the 
sweeping out of small pores during grain growth on recrystallization [181420] . 
The rim structure can be described taking into consideration many specific 
characteristics of the fuel [18J but in this paper we are not going to touch 
these interesting questions. 

Here we study the temperature distribution in the nuclear fuel rod of 
the reactor taking the influence of the high burn- up into account. Using the 
numerical method based on the difference equation of the heat conductivity 
we prove that the stationary behavior of the nuclear reactor is reached for the 
time less then one minute. This fact allows us to consider the temperature 
distribution using the stationary cases of the heat conductivity in the nuclear 
fuel reactor. At solution of the task we do not use the modern approaches for 
the nonintegrable differential equations (see [21H23] ). The task is solved by 
us taking the simple integration of equations. As a result we find the solution 
for the temperature distribution in the nuclear fuel rod in the analytical form 
for the stationary behavior of the reactor. Exact solutions allow us to analyze 
the temperature distribution in details and to evaluate the influence of the 
rim - layer and the zirconium oxide on the temperature distribution in the 
nuclear fuel rod. 

This paper is organized as follows. In section 2 we present the statement 
of the problem and we use the numerical methods to obtain the energy flow 
on the surface of cladding. Using the numerical simulation for solving task in 
section 3 we observe that the stationary case of the temperature distribution 
in the nuclear fuel rod arises very quickly in comparison with time of the 
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product of electricity production. So we can consider the stationary temper- 
ature distribution in the nuclear fuel rod. Using this fact in section 4 we give 
the stationary statement of the problem and in section 5 we find the exact 
solution of the temperature distribution in the nuclear fuel rod taking rim - 
layer and zirconium oxide into account for high burn-up fuel. In section 6 we 
discuss the particularities of the stationary thermal regimes in nuclear rods. 
As this takes place we consider the different cases of nuclear fuel rods: with 
gap and without gap; in the case of the burn-up and without rim - layer and 
zirconium and so on. We can allow to do it using the analytical solutions. 

2 The statement of the problem for the tem- 
perature distribution in nuclear fuel rod 
with rim — layer and zirconium oxide 

The geometry of the nuclear fuel rod is given in Fig. [TJ We take this geometry 
into account in this paper to study the temperature distribution in the nuclear 
fuel rod. 

Let us introduce the following notation: 2\ (r, t) is dependence of tem- 
perature on r and t in the fuel that is UO2, T2(r,t) is the dependence of 
temperature from radius and time in the gap, that is He, T 3 (r,t) is depen- 
dence of temperature on r and t in the cladding, T\ (r, t) is temperature in the 
rim layer, T 3 (r, t) is the temperature in the zirconium oxide, r is the radius 
of the central hole, Rq is the radius of the nuclear fuel, Ri is the external 
radius of the gap in the fuel rod, R2 is the radius of the cladding, Rq is the 
internal radius of the rim - layer, R2 is the external radius of the zirconium 
oxide. 

Nonlinear differential equation for the description of the temperature in 
the nuclear fuel rod can be written as 

dT 1 d ( dT\ 
C p (r,T) P(r,T)— = -— ^\(r,T)r—j +q(r), r < r < R 2 , (1) 

where C p (r, T) is the specific heat capacity, p(r, T) is the mass density, A(r, T) 
is the thermal conductivity, q(t) is the uniform volumetric heat generation 
rate which is independent of the temperature. 

We take the initial condition for the temperature in the form 

T(r,t = 0) = ip . (2) 

We used different initial conditions for the numerical simulation taking into 
account the initial behavior of nuclear reactor. 



3 




Figure 1: Geometry of the nuclear fuel rod: r is central hole, Rq is radius 
of the fuel, R is the radius of the rim layer, i? x is external radius of gap, R 2 
is radius of ZrC>2, R2 is external radius of Zr. 



The boundary values for the numerical simulation we took at r = r in 
the form 



dT 

Or 



= 0, (3) 

r=ro 



and at r = R 2 

T(r = R 2 ,t)=T w , (4) 

where T w is the temperature of the wall. 

Heat capacity (C p ) of oxide fuels is very important parameter for evalu- 
ation of fuel temperature at normal, transient and accidental conditions of 
light water reactor pQ. Heat capacities of undoped and impurity- doped UO2 
pellets have been measured by many researchers. The values for densities 
and heat capacities of nuclear fuel, rim layer, gas in gap and in cladding are 
different from each other. We take formulae for these values as sectionally - 
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continuous function 



p(r,T)C p ( ri T) 



, p l (T)C p>1 (T), ifr <r< J R , 

fc(T)C P)1 (T), ifi?o<r< J R , 

p 2 C P:2 , tfR <r<R ll 

p~ 3 C P:3 , if Ri < r < R 2 , 

p 3 C Pt3 , if R 2 < r < R 2 . 



We have taken physical parameters and dependencies of the densities and the 
specific heat capacities on temperature for the fuel, the gap and the cladding 
from the papers [HE1IZ] ■ 

Heat conductivities for the fuel and for the rim layer could be expressed 
by a classical phonon transport model pH427] 



Ai(T) 



A + cbu + BT 1 



(5) 



where A, B and c are constants and bu is the local burn-up in the nuclear 
fuel and in the rim - layer. 

For numerical simulation of the temperature distribution in the nuclear 
fuel and in the rim - layer We assume that the coefficient of the heat con- 
ductivity is expressed by the sectionally - continuous function 



A(T, r) 



A+BT' 
1 

A1+B1T' 

Ai, 
A 2 , 
A 2 , 



if r < r < R 
if Ro < r < Rq 
if #o < r < Ri 
iiR l <r < R 2 
HR 2 <r< R 2 , 



where A, B, A\, B\, \\, \\, X 2 , X 2 are constants. For numerical simulation 
of task we have taken values of these parameters from papers [2IH27]. 

The volumetric heat generation rate for the numerical simulation can be 
written as 

{<?o, if r < r < Ro 
9i, if R <r < R 
0, if Ro < r < Rq, 

where q and qi are constants. Task (JTj) - (J^J) is solved by us in the next 
section using the numerical methods. 
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3 The numerical simulation of the tempera- 
ture distribution in the nuclear fuel rod 

Consider the equation 

dT 1 d ( dT\ 
C p (r,T)p(r,T) — = --^X(r,T)r—j+q(r), r < r < R 2 , (6) 

Multiplying |6} on r and integrating with respect to r at constant q and 
qi, we have 

dE C R2 dT 

— = 2n J r\(r,T)^:dr + 7i q (R 2 -r 2 ) + n qi (R 2 -R 2 ), (7) 

where 

E = 2n C p (r,T)p(r,T)T(r,t)dr 

J r 

Let us denote the flow of energy via the surface unit as W, then Eq.([H]) 
can be written as 

BE 

— = -2-KRlW + >Kq{Rl- r 2 ) + n Ql (R 2 - R 2 ), (8) 

We obtain that the flow of energy for stationary behavior of the nuclear 
fuel rod takes the form 

w = Yr 2 ( q ~ r ° 2) + qi {r2 ° " ^° 2) ) ' (9) 

To solve the task ([T]) - (jlj) we use the numerical simulation. With this 
aim we introduce the grid on coordinate r and time t. We take points on r 
and t using the formulae 

rj = jh, (j = 0,...J); t n = rn, (n = 0,...N), (10) 

where h and r are steps on coordinate and time, J corresponds to the radius 
i?2, N corresponds to the stationary state of the nuclear fuel rod. 

Let us introduce the grid functions at t — t n and r — rj. In this case we 
have 

T(r v t n )-T^ qj~q(r 3 ,t n ), 

(11) 

a* ~ a(t;, rj ), c p (t;\ rj ) ~ c; d , p (t^ Tj ) ~ P ]. 
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We change the differential operators in Eq.([T]) on difference operators 
using the formulae 



dT 
~dt 



rpn+l rpn 

3 3 



r=hj;t=Tn+i 



dT 

Or 



rpn+l rpn+1 



r=hj+^;t=T(n+l) 



(12) 



Taking the grid functions ( II ip and the approximation of differential rela- 
tions (fl2"]l we have the following difference equation which is approximately 
equivalent to (JTJ 



rpn+l rpn 

3 3 



A 



Ti \ h 2 ) 



-A n 1 

3-2 



3-1 



71+1 



rpn+l 



h 2 



(13) 



(i = l,...J-l), 



where we use notation 

Difference equation (Tl3|) can be written as 

/in ™+l n" T"+l i Tjn rpn+l rpn 



(j = l,...J-l). (14) 



(15) 



where 



T 

V 2 



r 3+l 



2^ 



(16) 



d: 



T 
h 2 



3+1 



1 + ^j+1 

9 r 1 ™ n' 



r 

^2 



^7-1 



X] + X]^ 



o v k f>" 



(17) 



r 

7^ 



r.i-i 



(1* 



(19) 
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Difference equations (fT5|) are the system of the algebraic equations which 
corresponds to the differential equation ([I]). This system of equations allows 
us to find the temperature distribution taking the boundary values and initial 
conditions. 

We assume that the temperature distribution on time t n = nr is known. 
The task is to solve the system of the algebraic equations and to find the 
temperature distribution at the moment t n+1 = (n + l)r. 

From the initial condition we know that T(r, t — 0) = tpo and conse- 
quently we know the temperature distribution on the radius at n = 

if = ifo (20) 
Taking into consideration the boundary conditions ([3]) and (jlj), we have 
T n = T n^ Tj =Ti, (n = 0,...,N) (21) 

Using the temperature distribution on the first layer on time (1201) and 
conditions (1211) . we solve the system of the algebraic equations ffT5l) and find 
the temperature on the first layer at t l = r. Then taking into consideration 
values Tj, we obtain values of temperature T 2 at t 2 = 2r and so on. 

The system of algebraic equations f[T5"j) can be solved by the sweep method. 

On figure [2] we can see the evolution on temperature in nuclear fuel rod in 
the case qi = 200 W/cm. Here and later qi = q ir (i?Q — rfy+qi % (R^ — Rq), 
where q\ = 2q . We have used the initial condition for the temperature 
assuming that this temperature is equal to temperature of the wall 

T(r, t = 0) = T w (22) 

Other parameters of the mathematical model are used from the table 1. 

From figure [2] we can see that the stationary state of the nuclear fuel rod 
with burn-up is reached during 30 secund of time. 

Figure [3] demonstrates the output the temperature to the stationary 
regime of the nuclear fuel rod at the different values of the volumetrical 
source rate. We can observe that the time of the output on the stationary 
state is not more then 30 seconds. Therefore the basic temperature regimes 
of the nuclear fuel rods are the stationary regimes and we can consider the 
temperature distribution in the nuclear fuel rod as the stationary task. 

The numerical simulation of dependence for the heat flow on time at 
r = i? 2 is given on figure HI From this figure we can see that the time from 
the initial state to the stationary behavior is 40 seconds. This time depends 
on the value of source and the time grows with increasing source. However 
this time is not more 45 seconds for available sources in nuclear fuel rod. 
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Figure 2: Evolution of temperature in nuclear fuel rod for high burn-up for 
q L = 200 W/cm. 1: t = s; 2: t = 1 s; 3: i = 2 s; 4: i = 3 s; 5: t = 5 s; 6: 
t = 10 s; 7: t = 30 s; 

4 The statement of the stationary problem 
for the temperature distribution in nuclear 
fuel rods with the rim - layer and the zir- 
conium oxide 

Using the numerical simulation we have obtained in the previous section that 
the basic behavior of the nuclear fuel rood is stationary. Consequently we 
can consider the temperature distribution in the nuclear fuel rod taking the 
above mentioned problem when ^ = 0. In this case we have the following 
statement of the problem for the temperature distribution 

£G™§) + '-"- •»<'<* (23) 



i d 



dTt 



rdr \A 1 + B{T X dr 



+ q x = 0, Rq < r < Rq 



(24) 



2600 




r, mm 

Figure 3: Evolution of temperature in the nuclear fuel rod for high burn-up 
at q L = 400 W/cm. 1: t = s; 2: t = 1 s; 3: t = 2 s; 4: i = 3 s; 5: i = 5 s; 
6: t = 10 s; 7: i = 20 s; 8: t = 40 s; 



Id/, dTo . 

-~T rX ^~r = °> i?o < r < i?i (25) 



r dr \ dr 



i^frAa^k ) =0. /?, < / /? 2 (20) 



-T- f rA 2^ ) = 0, ^2 < r < i? 2 (27) 



r dr \ dr 

where T\ (r) is the dependence of the temperature on radius in the fuel, T\ (r) 
is the dependence of the temperature on radius in the rim layer, T 2 (r) is the 
dependence of the temperature on radius in the gap, T 3 (r) is the dependence 
of the temperature on radius in the zirconium oxide, T 3 (r) is the dependence 
of the temperature on radius in the cladding, q and q\ are volumetric sources 
in the fuel and in the rim - layer. 
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Figure 4: Dependence of the heat flow on time via surface at r = R 2 for 
sources 1: q L = 100 W/cm; 2: q L = 200 W/cm; 3: q L = 300 W/cm; 4: 
q L = 400 W/cm; 5: q L = 500 W/cm; (6m = 120) 



To obtain the solution of the system of equations fl23|) - (j2Zj) we have to 
add the boundary values for this system of equations. In this case we use 
the boundary conditions at r = r , r = R , r = R , r = Ri r = R 2 and at 
r = R 2 . 

The boundary value for the temperature T\(r) at r = r$ takes the form 
dTi 



dr 



0. 



(2f 



r=r 



We have conditions at r = Rn in the form 



A + STi dr 



r=iio 



Ax + 5iTi rfr 



r=R 



(29) 



T 1 (r = J R ) = T 1 (r = J R ). 
In the case r = R®, we use the following conditions 



1 



At + B x Ti dr 



Ax 



dr 



Tx(r = # ) = T 2 (r = i? ). (30) 



r=R 
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We take the conditions at r = Ri 
dT 2 



Ai 



dr 



~ dT 3 
A2 



r=Ri 



dr 



T 2 (r = R 1 ) = f 3 (r = R 1 ). 



(31) 



r=R 1 



Assuming the flow at r = R 2 we obtain the conditions 

dT 3 



A ^ 

2 



A, 



r=R 2 



dr 



T 3 (r = Ro) = T 3 (r = R 7 



(32) 



r=R 2 



We believe that temperature on the surface of the cladding at r = R 2 is 
constant and this temperature is equals to the temperature on the wall T w [6] 



T*(r = R 2 



T 



(33) 



System of equations (j23jl - (JzTjl and conditions fl28l) - fl33l) describe the 
stationary temperature in the nuclear fuel rod. 



Exact solutions for the temperature distri- 
bution in the stationary nuclear fuel rod 
with rim - layer and film of zirconium ox- 
ide 



Solutions of the system of equations fl23|) - - fl33|) can be obtained in the 
analytical form. Integrating equation (12"3"|) with respect to r, we have 



dTi qr 2 

1 — — — w, 



A + 5Ti dr 

where C\ is arbitrary constant. 
Eq.(JMD can be written as 



(34) 



d\n(A + BT 1 ) Bqr _BC X 
dr 2 r 

Integrating Eq.(l3"oT) with respect to r again we obtain 
In (A + BTi) = ln(C 2 ) + J BC 1 ln (r) - 



(35) 



(36) 
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where C% is arbitrary constant as well. From Eg. ( 1361) we have the solution of 
Eq.(J23J 

m / s Co R r~i -B qr 2 A ,„_, 

ri(r) = jr BCl e— --. (37) 

By analogue we find the solution of Eq. f[24l) for dependence temperature on 
radius in rim layer 

f,(r) = Jr**.^-^, (38) 

where Ci and C2 are constants of integration as well. 

The general solutions of Eqs. ( 1251) - fT2"T|) are determined by formulae 



T 2 (r) = ^ln(r) + C 4 , (39) 

T 3 (r) = Sln(r) + C 6 , (40) 
A2 

T 3 (r) = ^ln(r) + C 6 , (41) 
A2 

where (7 3 , (7 4 , (7 5 , (7 6 , (7 5 , and C 6 are arbitrary constants as well. 

In order to find the temperature distribution in the nuclear fuel rod we 
have to obtain values of the arbitrary constants Ci, C2, Ci, C2, C3, C4, C5, 
C*6, C5, and Cq taking conditions f l28|) - fl33l) into account. 

Taking the boundary value into account 



dTi 



dr 



0, (42) 



r=ro 

we obtain the constant C\ in the form 

C X = (43) 

Substituting solutions ( )3*7|) - fT4Tl) into conditions fl28|) - fj33|) . we have the 
system of algebraic equations with respect to constants (7 2 , Ci, (7 2 , C 3 , C 4 , 
(75,6*6,(75, (7 6 . This system of equations takes the form 



C2 5^ 
~B R ° 



BqR% 



C2 



.4 
5 



B1C1 



Kg 



^1 



(44) 
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Qro _ qRo = &i_ 9i^o ^ 



§L<^e-^-^ = ^lni? + C 4) (46) 



i? 2 i? 



^lni?! + C 4 = Slni?! + C6, (48) 
Ai A 2 



C 3 = C 5 , (49) 



Slni^ + Q; = ^ln£ 2 + C 6 , 



A 



A 2 



(50) 



C5 — C$, 



(51) 



A 2 



lni? 2 + C 6 = T U 



(52) 



Solving this system of equations with respect to constants C 2 , Ci, C* 2 , C3, 
C4, C5, C*6, C5, 6*6, we obtain dependencies of them on parameters of task 
in the form 



C 2 = B { T w + ^ + 




1 

*2 



^1 




B 1 g(r§-fl§)+-B 19l fl§ 



£1 



4 
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(53) 



Ci = 



g ~ g) | gi (g 



(54) 



C 2 = ^ BxT^ + Ai + BiCaln 



-R2 \ 2 (#A^ 



(55) 



? (rg - flg) gi (*8 - *g) 

C 3 = " A — : L + : ^ 



(56) 



C4 — r^, 



' g^g-flg) | gi(^-^ 



In 



(57) 



. 9 (r 2 - flg) 9l (flg - i2g) 
C 5 = ^ + — ^ — ^ 



(58) 



— T w 



In 



#2 



(59) 



C, = 



3 1 



"g ~ g) | gi (g ~ g) 



(60) 



C 6 



T — 



' g(r 2 -^) | gi(g-g 



In I R" 2 



(61) 
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Substituting these constants into the general solutions of the system of 
equations we have formulae to determine the temperature distribution in the 
nuclear fuel rod. The temperature distribution in the fuel without rim - layer 
takes the form 



T 1 (r)= {T w + ^- + W ln 



r \ A i / Ri \ *2 / R 



Ri 



Ro 



Ro 



Ro 

Rq 



B 1 W 1 iBr, 

r 



Rq 



2 q^^Rl-R^ + qBjRl-r 2 ) 



A A 1 



2 qB(Rl-r 2 ) A 
P 4 



where we use the notation 



Q (r 2 o - Rl) qi (Ri - Ri) q U ~ %) 



(62) 



+ 



qiRj (63) 



The temperature distribution in the rim - layer has dependence on pa- 
rameters of task and on the radius in the form 



Ti(r) = {T w + ^- + W \n 



r \ A i 

Ri 



Rl \ A 2 / i?2 



Ro 



1 1 

*2 



Ro 



Rq J Bi 

The temperature distribution in the gap can be written as 



(64) 



T 2 (r)=T w + W \n 



r \ x i 

Ri 



R\ \ A 2 / R 2 



Ro 



1 

A 2 



R2 



(65) 



The temperature in the film of the zirconium oxide is expressed by formula 



f 3 (r)=T w + W In 



r \>2 / Ro 



Ro 



Ro 



(66) 



The temperature distribution in the cladding is found by formula 



T 3 (r) = T w + W In 



T \ A 2 

R~2 



(67) 
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The temperature distribution in the nuclear fuel rod is described by the 
continuous function with the break of the first derivative 



T(r) 



2i(r 
fx(r 
T 2 (r 
T 3 (r 
T 3 (r 



if r < r < Rq, 
if Rq < r < Rq, 
if Ro < r < Ri, 
if R x < r < R 2 , 
HR 2 <r < R 2 . 



(68) 



Using formula ( |68|) we can find the temperature distribution in the nuclear 
fuel rod at given volumetric source rate in the fuel q and in the rim - layer q\. 
We can find the temperature distribution in the case of different thickness of 
the film for the zirconium oxide and at different sizes of the rim - layer. From 
formulae f )62|) - (j67|) we can obtain a number of partial dependencies for the 
temperature distribution in nuclear fuel rod. We can find the temperature 
distribution in the nuclear fuel rod without rim - layer, without the zirconium 
oxide, without gap and so on. With this aim we have to take R = Ro (rim 
- layer is absent), R 2 = R 2 (the zirconium oxide is absent), R = Ri (gap is 
absent) and so on. 



6 Results and discussion 

Let us consider the temperature distribution in a nuclear fuel rod taking the 
rim - layer and the film of zirconium oxide into account. To analyze the 
temperature distribution we use sizes of nuclear fuel rod from the table [TJ 



Table 1: Parameters of the nuclear fuel element 



ro, mm 


R , mm 


Ri, mm 


R 2 , mm 


0.75 


3.775 


3.865 


4.550 



Parameters of heat conductivities of the fuel, the helium in the gap, the 
film of the zirconium oxide and the cladding of the zirconium are given in the 
table [2J We assume, that the heat conductivity is 2 times less in a rim — layer, 
than the heat conductivity in the fuel core. 



Table 2: Parameters of heat conductivities of fuel element 



A, mm-K/W 


B, mm/W 


Ai, W/mm-K 


A 2 , W/mm-K 


A 2 , W/mm-K 


T K 


43.8 


0.2294 


0.3 • 10" 3 


1.8 ■ 10" 3 


22.0 • 10~ 3 


620 
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We assume that the temperature of the coolant T w is equal to 620 K 
350 °C). 
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Figure 5: Dependence of temperature on radius at qi = 200 W/cm for cases: 
1: gap and cladding are out; 2: gap is absent; 3: bu = 0; 4: few = 120; 

Some thermal regimes of temperature distributions in the nuclear fuel 
rod are given on figure [5] and on figure [6] by the formulae (162]) — (|67|) . The 
first figure illustrates curves for the value of the source qi = 200 W/cm. 
The second figure demonstrates the thermal regimes in the case of the source 
q L = 400 W/cm. 

Curve 1 illustrates the temperature distribution when the gap and the 
cladding are absent. This case is possible when there is breakage of the 
cladding and the coolant has contact with the fuel. Temperature of the fuel 
in this case is less then in other cases. The decreasing of temperature derives 
from the fact that energy of the fuel goes to the coolant because the heat 
conductivity of fuel is more then the heat conductivity of the helium in the 
gap. 

Curve 2 in figures corresponds to the case with cladding but when the gap 
is absent. One can observe that the temperature is higher then in previous 
case but this one is less then the temperature in the nuclear fuel rod for the 
normal behavior. 

Curve 3 illustrates the case of the normal behavior of the nuclear fuel rod. 
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Figure 6: Dependence of temperature on radius at source q^ = 400 W/cm 
for cases: 1: gap and cladding are absent; 2: gap is absent; 3: bu = 0; 4: 
bu = 120; 

The temperature distribution is higher of the previous cases but this one is 
less then in the case with the formation of the rim - layer. 

Curve 4 on figure |5] and on figure [6] corresponds to the thermal regime in 
nuclear fuel rod of high burn-up. We take into account the burn-up equal 
120. We can see that the temperature in this case of the nuclear fuel rod is 
much more then other cases. The comparison of the temperature in this case 
with the point of the melting fuel shows that there is dangerous behavior of 
the nuclear reactor because it is possible the melting of the fuel at formation 
of the rim - layer. This fact illustrates that in the case of increasing of the 
rim - layer we have to decrease the value of the source. 

At numerical calculations we suppose that the film of the zirconium oxide 
is formed at the expense of the gap and the cladding at various values of burn- 
up of the fuel. As this takes place the thickness of the gap is decreased. We 
assume that the film of the zirconium oxide takes 90 % of the cladding and 
10 % of the volume in the gap. 

The dependence of temperature on radius in the case of the normal be- 
havior in nuclear fuel rod and for high burn-up is given in figure [7] (q = 0.2) 
and in figure E] (q=0.4). 
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Table 3: Parameters at burn- up of fuel 



Burn-up, GWt • twenty four hours/t-U 





60 


80 


100 


120 


Thickness of a layer Zr0 2 , mm 





0.0168 


0.0224 


0.028 


0.0336 


Thickness of a rim — layer, mm 





0.1 


0.15 


0.2 


0.25 



1500 
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Figure 7: Dependence of temperature on radius at qi = 200 W/cm for cases: 
1: bu = 0; 2: bu = 120 (q L = 200 W/cm); 3 : bu = 120 (q L = 227 W/cm) 

We can observe that the high burn-up leads to the increasing of the 
temperature in nuclear fuel rod. The point of the melting of UO2 is equal to 
3000 K. In the case of the high burn- up equal to 120 we can believe that the 
fuel can melt and we need to decrease the power of the source in the nuclear 
fuel rod. 

Dependence of temperature on source q in the point r = r of the nuclear 
fuel rod is presented in figure M at various depths of burn-up of the fuel. 

From this figure we can see how should we have to increase the power 
in nuclear fuel rod to have the normal behavior of the fuel element. The 
formation of the rim - layer leads to increasing of the heat conductivity. As 
this fact takes place the power of the source in nuclear fuel rod is decreased. 
The linear power of the nuclear fuel rod is decreased as well. Temperature 
in the fuel and in the rim - layer is decreased. We can see this process in 
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Figure 8: Dependence of temperature on radius at q — 300 W/cm for cases: 
1: bu = 0; 2: bu = 120 (g L = 4 W/cm); 3: 6m = 120 (q L = 453 W/cm) 

figure U\ and in figure [HJ Therefore we have to increase the power in nuclear 
fuel rod. 

7 Conclusion 

Let us shortly formulate results of this paper. We have studied the temper- 
ature distribution in nuclear fuel rod taking into account the fuel, the rim - 
layer, the gap, the film of the zirconium oxide and the cladding. Solving the 
nonstationary task of the temperature distribution in the nuclear fuel rod by 
the numerical simulation we have obtained that the evolution of the temper- 
ature in the rod to the stationary behavior goes the short time. This time is 
less then 45 seconds for the maximum power in nuclear fuel rod. This fact 
allowed us to consider the solution of the stationary behavior of the nuclear 
reactor. We have solved this task using the analytical method and we have 
found the exact solution of the temperature distribution in the nuclear fuel 
rod. We analyzed the different thermal regimes for the stationary behavior 
reactor and have shown the important role of the rim - layer of high burn-up 
nuclear fuel rod. 
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Figure 9: Dependence of temperature in point r = r on linear power of 
nuclear fuel rod (q\ = 2q) at different burn-up 1: bu = 0; 2: bu = 60; 3 : 
bu = 120; 
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